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Abstract. In this article we try to bridge the gap between the quantum 
dynamical semigroup and Wigner function approaches to quantum open sys- 
tems. In particular we study stationary states and the long time asymptotics 
for the quantum Fokker-Planck equation. Our new results apply to open 
quantum systems in a harmonic confinement potential, perturbed by a (large) 
sub-quadratic term. 



1. Quantum Fokker-Planck model 

This paper is concerned with the mathematical analysis of quantum Fokker- 
Planck (QFP) models, a special type of open quantum systems that models the 
quantum mechanical charge-transport including diffusive effects, as needed, e.g., in 
the description of quantum Brownian motion [T] , quantum optics [5] , and semicon- 
ductor device simulations [3]. We shall consider two equivalent descriptions, the 
Wigner function formalism and the density matrix formalism. 

In the quantum kinetic Wigner picture a quantum state is described by the real 
valued Wigner function w(x, v,t), where (x,v) G R 2 denotes the position- velocity 
phase space. Its time evolution in a harmonic confinement potential Vo(x) = 
with uj > is given by the Wigner Fokker-Planck [5J Q] [5] (WFP) equation 

(1.1) dtw = u> 2 xd v w — vd x w + Qw , 

Qw = 2 r yd v (vw) + D pp A v w + D qq A x w + 2D pq d v d x w . 

The (real vauled) diffusion constants D and the friction 7 satisfy the Lindblad 
condition 

(1.2) A:=D pp D qq -Dl q - 1 2 /A>Q, 

and D pp , D qq > 0. Moreover we assume that the particle mass and h are scaled to 
1. 

WFP can be considered as a quantum mechanical generalization of the usual 
kinetic Fokker-Planck equation (or Kramer's equation), to which it is known to 
converge in the classical limit h — ► 0, after an appropriate rescaling of the appearing 
physical parameters [8] . The WFP equation has been partly derived in Ref . [9] 
as a rigorous scaling limit for a system of particles interacting with a heat bath of 
phonons. 

In recent years, mathematical studies of WFP type equations mainly focused 
on the Cauchy problem (with or without self-consistent Poisson-coupling) [THl HH 
EJ [13j [12l [14] . In the present work we shall be concerned with the steady state 
problem for the WFP equation and the large-time convergence to such steady states. 
Stationary equations for quantum systems, based on the Wigner formalism, seem 
to be rather difficult. For a purely quadratic confinement potential, this problem 
was dealt with in Ref. [5] using PDE-tools. The extension to harmonic potentials 
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with a small, smooth perturbation was recently obtained in Ref. [T5] using fixed 
point arguments and spectral theory. Here we consider large perturbations of the 
harmonic potential. To this end we shall work in the density matrix formalism, 
using tools from operator theory. 

In the density matrix formalism a quantum state is described by a density matrix 
p € the cone of positive trace class operators on some Hilbert space h. Its 

time evolution is governed by the linear QFP equation or master equation 

(1-3) ^ = 

with the Lindbladian 

= ~\ [p 2 +LU 2 q 2 + V(q),p] -ij[q,{p,p}] 

(1-4) - Dgg[p,\p,p\}-D pp [q,[q,p}}+ 2D pq [q, [p, p}} , 

where V(q) is the perturbation of the harmonic potential. 

Global in time solutions to such master equations were established in Ref. [TS| 
(nonlinear QFP-Poisson equation) starting from the construction of the associated 
minimal quantum dynamical semigroup (QDS) [16j . 

General methods for the study of quantum master equations and their large time 
behavior, including the existence of steady states and convergence towards them 
were developed in Quantum Probability. 

Applicable sufficient conditions for proving uniqueness, i.e. trace preservation, 
of the solution obtained by the minimal semigroup method were given in Ref. [17j 
(see also Ref. [H]). A criterion based on a non-commutative generalization of 
Liapounov functions for proving the existence of steady states was developed in 
Ref. [H]. The support of steady states and decomposition of a quantum Markov 
semigroup into its transient and recurrent components were studied in Ref. |20j and 
[21]. When the support of a steady state is full, i.e. it is faithful, uniqueness of 
steady states and convergence towards steady states can be deduced from simple 
algebraic conditions based on commutators of operators appearing in a Lindblad 
form representation of the master equation (see [22] for bounded and [53] for un- 
bounded operators). Many of these methods generalize those of stochastic analysis 
in the study of classical Markov semigroups and processes. We refer to the lecture 
note [U] for a comprehensive account. 

In this paper we study the master equation (|1.3p by the above methods. We first 
prove the existence and trace preservation (i.e. uniqueness) of solutions and then 
the existence of a steady state. If the diffusion constants D pp , D pq , D qq , and the 
friction 7 satisfy the Lindblad condition (|1.2j) with the strict inequality, we prove 
that this quantum Markovian evolution is irreducible in the sense of Ref. [25] . As a 
consequence, steady states must be faithful and one can apply simple commutator 
conditions on the operators in the GKSL representation to establish uniqueness of 
the steady state and large time convergence towards this state. 

When A = we conjecture that (see Sect. H]), unless V is zero and the limiting 
conditions D qp = —^D qq , D pp = uj 2 D qq are satisfied, the quantum Markov semi- 
group is still irreducible. But the invariant subspace problem that has to be solved 
for proving this becomes very difficult and we were not been able to solve it. 

The paper is organized as follows: In Section [2] we review the equivalence of 
the kinetic Wigner formalism and the Lindblad approach to open quantum sys- 
tems. Some technical preliminaries are presented in S}3] and [H In S}5] we construct 
the minimal QDS for (|1.3p . (|1.4[) with external potentials that grow at most sub- 
quadratically. The markovianity of the semigroup is proved in This yields 
uniqueness and mass-conservation of the solution to (|1.3|) . (|1.4p . In |J7]we establish 
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the existence of a steady state and in $3] we prove that the solution converges to this 
unique steady state for arbitrary initial data provided the Lindblad condition (|1.2p 
is fulfilled with strict inequality. The limiting case A = is studied in Sj9] 



2. Passage from the Wigner equation to the master equation 

In this section we show how to pass from the Wigner language to the GKLS 
(Gorini, Kossakowski, Sudarshan [55]; Lindblad [57]) language. In order to keep 
the notation simple, we shall confine our presentation to the one dimensional case. 
However, the results extend to higher dimensions. The underlying Hilbert space of 
our considerations is h = L 2 (R) . We denote by q and p the standard position and 
momentum operators (p = —id x ). They satisfy the canonical commutation relation 
(CCR) [q,p] =il. 

The Wigner function w(x,v,t) of a state pt = % t {p), is (up to normalization) 
the anti Fourier transform of 

(2.1) ^(e,r ? )=tr( Pt e- i ^+^) . 

Using (|2.ip we shall now transform the WFP equation (jTTTJ) into an evolution 
equation for the corresponding density matrix p. 
As a consequence of the CCR we have 

(2.2) e -H.£q+vp) = e -»£9 e -iw e »£»7/2 ; 

(2.3) e -i(iQ+vp) — e -ivv c -i£,i c -i(,n/^ ; 

(2.4) Q-iiii+vv) — g-i»7P/ 2 e -*£? e -"7P/ 2 . 



Assuming that p is sufficiently regular, by differentiating ()2.2p and (|2.3I) and using 
the cyclic property of the trace, we find 

%<p(t,V) = -*tr(pge-^ + "^) + % -£<pfov), 
9MZ,V) = -itr(gpe-^+^)-|^,r 7 ). 
Subtracting, and respectively, summing the above equations we have 



In a similar way we obtain formulae for the products and derivatives with respect 
to £: 

&(Z,V) = tr{[p,p]e-^+^) , 
fyvfov) = -^tr({p,p}e- i ^ + '^) . 

The Wigner function is the anti Fourier transform of tp 

w(x,v) = (J^j J^S^+^cp^^dr,. 
The factor is chosen such that the total mass is given by 

m = tr (p) = ip(Q, 0) = / w(x,v)dxdv . 
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Integrating by parts (and again assuming sufficient regularity and decay) we obtain 
1 



xw(x, v) 



2tt 



K 2 



i / -ifd^'+wAtpforiWdri 



2tt 



i 

2^ 



In a similar way one can calculate 

27T . B 2 \ 2 



d£dr/ . 



vw(x, v) = / e ^ x+T > v > ( itr ( {p, p) c -^ q+r ^ ) ) d£dr) , 

271" JR2 \ 2 V 7 / 

d x w(x,v) = ^- y e i ^+"^ (tr(i[p,p]e-^« 9+ " p) )) d^dry, 

a„w(x,tO = —[ e l ^+^) (tr(-i[q,p]e~^ q+r > p A) d£,dr) . 

2n J R 2 V V / / 

The above formulae lead to the following dictionary for translating a master equa- 
tion from the Wigner function language to the GKSL language: 



transformation on w 


transformation on p 


xw 
vw 
d x w 
d v w 


H?, p} 
\{v,p\ 

i[PiP\ 

-i[q,p] 



For the terms appearing in the WFP equation we have: 



transformation on w 


transformation on p 


xd v w 


-f{e, k,p]} = -f[s 2 ,p] 


vd x w 


\{p,\p,p}}=\\p 2 ,p] 


d v (vw) 


-|[?,{P.P}] = -|{P, [M) + P 


A v w 


-[?> hp}} 


A x w 


-\p, hp}} 


d x d v w 


\p, [q,p\] 



Using this dictionary we find the following GKSL form = £*(pt)) of the 
linear QFP equation 

(2.5) C,(p) = -\ [p 2 +u 2 q 2 ,p] -il[q,{P,p}} 

- D qq [p, [p, p]] - D pp [q, [q, p\] + 2D pq [q, [p, p]] . 

This corresponds to choosing A = p = 7 in (3.8) of Ref. [5] (see also [2E])- 
The dual equation of (|2.5p with an added perturbation potential V reads 



C(A) = l - [p 2 +u 2 q 2 + 2V(q),A] +i 1 {p,[q,A}} 

D qq [p, [p,A]] - D pp [q, [q,A]} + 2D pq [q, [p,A]] , A e #(h) . 
It can be written [I] in (generalised) GKSL form like 

1 2 

(2.6) C(A) =i[H,A]--^2 (L* e L e A - 2L\AL t + AL\L t ) 

1=1 
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with the "adjusted" Hamiltonian 

H = - {p 2 + uj 2 q 2 + j(pq + qp)) + V{q) 
and the Lindblad operators L\ and L2 given by 



(2.7) L x = p+y/2D pp q, L 2 = = p. 

V ^^vv V pp 

3. Key inequalities for the existence of a steady state 

We aim at applying the criterion for existence of a normal invariant state [19j 
by Fagnola and Rebolledo. To this end we have to find a positive operator X and 
an operator Y bounded from below, satisfying 

C(X) < -Y, 

which in addition both have finite dimensional spectral projections associated with 
intervals ] — 00, A]. To illustrate the technique we first present the computation for 
the harmonic potential only. The perturbation potential V(q) will be added later 
on. 

Consider the Lindbladian 

£*{p) = -o [p 2 + ^V,p] - i-f[q, {p,p}} 



2 

D qq [p, [p, p}} - D pp [q, [q, p}} + 2D pq [q, [p, p]] 



with dual [28] 



C(X) = l -[p 2 +L0 2 q 2 ,X]+t 1 {p,[q,X]} 

- D qq [p, [p, X]] - D pp [q, [q, X]] + 2D pq [q, ]p, X]] . 
Straightforward computations with the CCR [q,p] — it yield 
Lemma 3.1. The following formulae hold for f , g smooth: 

£(f(p)) = - t ^(qf'{p) + f'{p)q)-2 1 pf'{p) + D pp f"{p) 1 

C{g{q)) = \(pg'(q)+9'(q)p) + D qq g"(q), 

C(pq + qp) = 2 (p 2 - u: 2 q 2 ) - 2 1 {pq + qp) + AD pq . 

This suggests looking for X, Y given by second order polynomials in p and q 
{i.e. f(p) = p 2 , g(q) = q 2 ). Therefore we start studying some algebraic properties 
of these operators: 

Lemma 3.2. For all r, s > such that rs > 1 the operators 

rp 2 — (pq + qp) + sq 2 , rp 2 + (pq + qp) + sq 2 

are strictly positive and have discrete spectrum. Moreover all spectral projections 
associated with bounded intervals are finite dimensional. 

Proof. Let ro, so > be such that ro < r, sq < s and roSo = 1- Then 

WroP - y/soq\ 2 = r p 2 - (pq + qp) + s q 2 > . 

It follows that 

rp 2 - (pq + qp) + sq 2 = \^/r^p - *Js^q\ 2 + (r - r )p 2 + (s - s )q 2 . 

Therefore the resolvent of rp 2 — (pq + qp) + sq 2 is dominated by the resolvent of a 
multiple (indeed min{ (r — ro), (s — So) }) of the number operator \(p 2 + q 2 — 1). 
Since the latter is compact, also the resolvent of rp 2 — (pq + qp) + sq 2 is compact 
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and self-adjoint. Hence it has a discrete spectrum that might only accumulate at 
0. 

The proof for the second operator is the same. □ 

We choose X of the form 
(3.1) X = rp 2 + (pq + qp) + sq 2 

and compute 

C{rp 2 + (pq + qp) + sq 2 ) = -2(2 7 r - l)p 2 - 2ui 2 q 2 

+(s - 2 7 - uj 2 r)(pq + qp) + 2rD pp + AD pq + 2sD qq . 

The required conditions on X and Y (i.e. X > 0, C(X) < —Y) hold if 

rs > 1 , 

Aid 2 (2jr - 1) > \s - 27 - uj 2 r\ 2 . 

Letting r and s go to infinity with s — 2j — to 2 r constant (that can be 0, for 
simplicity), it is clear that, when 7 > 0, we can find r and s large enough satisfying 
the above condition. We take, e.g. any r > (2-f)~ 1 and s = 2-f + ui 2 r since 

,2 



rs = 2-yr + uj l r l > 1 + — - > 1 
47 2 



Y will be chosen later in Theorem [21 



We now add the perturbation potential. Let V : M — > M. be a smooth function 
satisfying a growth condition like 

(3.2) \V'(x)\<g v (l + \x\ 2 ) a/2 

with gv > and < a < 1. Hence, this perturbation potential is strictly sub- 
quadratic. It gives rise to one additional term in C(X), namely: 

i[V(q),p 2 } = -(pV , (q)+V'(q)p) , 

i[V(q) 1 pq + qp] = -2qV (q) . 

Therefore we find now 

(3.3) C(rp 2 + (pq + qp) + sq 2 ) 

= -2(2jr - l)p 2 - 2u?q 2 + (s - 2-f - uj 2 r)(pq + qp) 



-r (pV'(q) + V'(q)p) - 2qV'(q) + 2rD„ 
Note that (due to the positivity of \e lt2 p ± e~ 1,2 V'(q)^ 



AD 



2sD, 



qq ■ 



+ \{V'(q)) 2 ^ < pV'(q) + V'(q)p < ep 2 + \(V (q)f 



Therefore, playing on the s and the bound on the derivative of V, we can find the 
needed inequality C(X) < —Y. This will be used in |J7]to prove the existence of a 
steady state. 



4. Domain problems 
First we define the number operator N 

Dom(iV) = I u G h Nu eh| = |ueh p 2 u 1 q 2 u 6 h J , 



|(p 2 + q 2 — 1) on h with 



where the last equality follows easily from ||iVu|| 2 < 00 by an integration by parts. 
C7?°(R) is a core for N (cf. Ref. [30], e.g.). Let X be the self-adjoint extension of 
(|3.ip (still denoted by X). Dom(X) is its maximum domain and C£°(R) is a core 
for X. 
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The position and momentum operators are defined on Dom(A f1 / 2 ). Both q and 
p have, by Nelson's analytic vector theorem, self-adjoint extensions that will be still 
denoted by q and p. 

First we shall compare the domains of ./V and X. To this end we need 
Lemma 4.1. Let r, s > with rs > 1 and define 

R := r l ' 2 p + r-Wq, S := s 4 ' 2 q + s^p . 
Then, for all u G (K) the following identities hold 



(u,X 2 u) = (s- ^ 2 (u,q L u) 



2 Is - i J (m, gi? 2 g u) 
( U;J R 4 W )-2(r S -l)|| U || 2 , 



(u,X 2 u) = - 2 (u,p 4 u) + 2 - {u,pS 2 pu) 

+ (u,5 4 u)-2(rs-l)||u|| 2 . 

Proof. Since u belongs to the domain of any monomial in p and q the proof can be 
reduced to the algebraic computation avoiding writing it's. 

Starting from the identity X = 9q 2 + R 2 with 6 = s — 1/rwe have 

X 2 = 9 2 q 4 + 9 (q 2 R 2 + R 2 q 2 ) + R 4 . 

The mixed product term can be written in the form 

q 2 R 2 + R 2 q 2 = qR 2 q + q[q,R 2 ] +qR 2 q+[R 2 ,q]q 

= 2qR 2 q + qR [q, R] + q [q, R] R + [R, q] Rq + R [R, q] q 
= 2qR 2 q + 2ir 1/2 qR - 2ir 1/2 Rq 
= 2qR 2 q + 2ir 1/2 [q, R] 
= 2qR 2 q - 2r . 

The conclusion is now immediate. □ 

The following lemma gives similar inequalities for the operators rp 2 + sq 2 (i.e. 
X without mixed products) that will be useful in the sequel 

Lemma 4.2. For all r, s > and u G C^°(R) we have 

(u,(rp 2 + sq 2 ) 2 uj > (r A s) 2 (u, (p 2 + q 2 ) 2 uj - 2 (rs - (r A s) 2 ) ||uj| 2 , 
(u, (rp 2 + sq 2 ) 2 uj < (r\/s) 2 (u, (p 2 + q 2 ) 2 uj + 2 ((r\/s) 2 -rs) \\u\\ 2 , 

where r A s = min { r, s } and rVs = max { r, s } 
Proof. Indeed 

/ 2 i 2\ 2 24, / 2 2 , 2 2\ , 24 

(rp + sq I = rp+rs(pq+qp)+sq 
= r p + 2rspq p + s q — 2rs 
> (r A s) 2 (p 4 + 2pq 2 p + q 4 ) - 2rs 

= (r A s) 2 (p 2 + q 2 ) 2 - 2 (rs - (r A s) 2 ) . 
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Moreover 



/ 2 i 2\ 2 24, / 2 2 , 2 2\, 24 

(rp + sg J = rp+rs[pq+qp)+sq 



= r 2 p 4 + 2rspq 2 p + s 2 q i — 2rs 
< (r V s) 2 (p 4 + 2pg 2 p + g 4 ) - 2rs 

= (r V s) 2 (p 2 + g 2 ) 2 + 2 ((r V s) 2 - rs) . 

This completes the proof. □ 

Proposition 4.3. The domains of the operators N and X coincide. 

Proof. We first show that Dom(X)CDom(iV). By Lemma T4.ll for all u £ we 
have 

\\Xu\\ 2 + 2(rs - l)j|u|| 2 > i min { (s - r" 1 ) 2 , (r - s" 1 ) 2 } (u, (p 4 + q 4 )u) 

> \ mm { (s r- 1 ) 2 , (r - s' 1 ) 2 } (u, (p 2 + q 2 fu) , 

where we used the elementary inequality (p 2 +q 2 ) 2 < 2(p 4 +q 4 ) and r, s > 0; rs > 1. 
Therefore we find a constant c±(r, s) > such that 

(4.1) ||X U || 2 + 2(rs - l)j|u|| 2 > Cl (r,s) \\(p 2 + q 2 )u\\ 2 . 

Now, if u € Dom(X), then there exists a sequence {u n ) n >i is converging in 
norm to u such that (Xu n ) n >i converges in norm to Xu (since X is closed). The 
above inequality shows then that the sequence (Nu n ) n >i is also norm convergent 
and its limit (N is closed) is Nu. 

This shows that Dom(X) C Dom(iV) and the two domains coincide since the 
opposite inclusion holds true by the construction of X itself. □ 

For future reference we briefly recall the definition of the annihilation operator 
a := -j=(q + ip) and the creation operator := ^(<Z ~ ip)- Using the isomor- 
phic identification of h with / 2 (No) (via the eigenfunctions of N, i.e. the Hcrmitc 
functions), these operators can also be represented as 



ae = 0, ae j+ i = y/j + lej, a^Cj = \/ j + le j+1 , for j£N . 
Here, (ej)j>o is the canonical orthonormal basis of Z 2 (No). 

5. Construction of the minimal quantum dynamical semigroup 

In this section we shall establish the existence of the minimal quantum dynamical 
semigroup (QDS) for the Lindbladian (|2.6p . To this end we first consider the 
operator G, defined on Dom(iV) by 

1 

(54) + (D pq - % -fj(pq + gp) + ^-iV(q). 

We suppose that the potential V is twice diffcrcntiable and satisfies the growth 
condition 



G = --(L* 1 L 1 +L* 2 L 2 )-iH = -(D qq + -)p 2 -[D pp + —)q 2 



Proposition 5.1. The domains of the operators N, G, and G* coincide. 

Proof. Since Dom(iV) C Dom(G) by construction, it suffices to check the opposite 
inclusion. To this end we proceed as before finding an estimate of ||-/Vit|| by the 
graph norm of G. Putting 

Go = -- (L{Li + L*L 2 ) = -D qqP 2 +D pq {pq + qp) - D pp q 2 + 1 , 
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we have for all u G G C °°(R) 

||G U || 2 = ||G QU || 2 + {u,i[H, G ] u) + \\Hu\\ 2 . 

A straightforward computation yields 

i[G ,H) = -2( 1 D qq + D pq )p 2 + 2( 1 D pp +u 2 D m )q 2 +2D pq qV'(q) 

+ (D pp -U 2 D qq ) (pq + qp) - Dqq (pV'(q) + V'(q)p) . 

The above commutator has quadratic monomials in p and q and other terms in- 
cluding V'(g) whose growth is sublinear due to our hypothesis on the potential. It 
follows then that we can find a constant C3 > such that 

\(u,i[G ,H]u)\ <c 3 (u,(p 2 + q 2 )u) . 

Hence, for all e > 0, by the Schwarz inequality 

|(u, i [Go, H] u)\ < e || (p 2 + q 2 )u\\ 2 + c 2 e- l \\u\\ 2 . 

It follows that 

(5.2) ||Gu|| 2 > !|G 0U || 2 - e ||(p 2 + g 2 )u|| 2 - c^lM) 2 . 

Case 1: 

Let D pq ^ and D pp D qq > D pq (cp. to the Lindblad condition l|1.2|) ). Then 

99 Dpp -Dpp^qq 

n n = n2 > 1 

^pq ^pq pq 



and 



( |% 2 - (pq + qp) + ^q 2 j + \ 



Therefore Go is a multiple of an operator like those of Lemma 13.21 (up to a con- 
stant). By inequality (|4.1|) we can find positive constants C4, C5 > such that 
||Goit|| 2 > C4 \\(p 2 + 9 2 )it|| — C5 ||u|| 2 . Therefore we derive the inequality 

(5.3) ||G U || 2 > ( C 4-e)||b 2 + g 2 ) U || 2 - C6 ||u|| 2 

for all u € C£°(R). Choosing e = C4/2, we find an inequality allowing us to repeat 
the above argument for a sequence (u n ) n >i (see the proof of Prop. FO)) and prove 
the inclusion Dom(G) C Dom(iV). 

Case 2: 



Let D pq = 0. The estimate (|5.3[) now follows from (|5.2|) and Lemma [4.21 We con- 
clude again that Dom(G) C Dom(TV). 

Case 3: 

Let D pq and D pp D qq — D pq (and hence 7 = 0). To recover ()5.3|) we start from 

(5.4) ||Gw|| 2 > \\Hu\\ 2 - e\\(p 2 + q 2 )u\\ 2 - cle-'Wuf 

(in analogy to (|5.2|) ). where we have now H = ^(p 2 +uj 2 q 2 ) + V(q). From Lemma 
we obtain 

\\(p 2 +q 2 )uf< Cl \\^{p 2 +u; 2 q 2 )u\\ 2 + C 2 \\u\\ 2 

<2 Cl (\\Hu\\ 2 + \\V(q)u\\ 2 ) + c 2 \\u\\ 2 . 
Since V(q) is sub-quadratic we have 

\\V(q)u\\ 2 < 4 A|| 2 + cje-'Wuf < e\\ (p 2 + q 2 )u\\ 2 + Ci \\u\\ 2 . 

Hence, 

(5.5) ||(p 2 + g 2 H| 2 -c 5 || U || 2 < C6 ||i/ M || 2 , 
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and we conclude by combining (|5.4p and (|5.5p . □ 

Remark 5.2. From Propositions 14.31 and 15.11 we infer that the graph norms of 
G, X and N are equivalent. Hence G is relatively bounded by X. 

The operator G is clearly dissipative because 

i 2 

R(u,Gu) = --^2\\L e u\\ 2 < 
£=i 

for all u £ Dom(iV). Therefore, by Prop. 3.1.15 of Ref. [29] it is closable and its 
closure is dissipative. We denote by the same symbol G the closure. Analogously, 
G* is also dissipative. 

Hence, the Lumer-Phillips theorem (see Theorem 3.1.16 of [23], e.g.) yields: 

Proposition 5.3. The operator G generates a strongly continuous contraction 
semigroup (Pt)t>o on h. 

Since Dom(G) = Dom(iV) c Dom(L^), t = 1, 2, and since 

2 

(Gv, u) u, v £ Dom(iV) , 

i=\ 

condition (H) of Ref. [Tjjj holds and we can construct T, the minimal QDS [H] 
associated with G and the L('s. 

6. Markovianity of the Quantum Dynamical Semigroup 

The hypotheses for constructing the minimal quantum dynamical semigroup 
with form generator 

2 

(6.1) £(A)[v,u] = {Gv,Au)+Y J ( L ev,AL e u) + {v,AGu), A G £)(h) 

l=i 

hold by Prop. PI 

Now we want to show that the minimal semigroup T is Markov and hence mass 
conserving. To this end we apply Theorem 4.4 in Ref. [17] , 

The algebraic computations in Sec. [3] suggest to consider an operator X of the 
form ([3~T]) with r > (27)" 1 and s = 2j + ui 2 r on the linear manifold C C °°(K). The 
algebraic computations of this section can be made in the quadratic form sense. 

We now check that the operator X satisfies the fundamental hypothesis C of 
Ref. |17j starting from domain properties: 

Proposition 6.1. The following properties of G, Li, and N hold: 

(1) Dom(G) = Dom(X) C Dom(X 1 / 2 ) and Dom(G) is a core for X 1 / 2 , 

(2) L £ (Dom(G 2 )) C Dompf 1 / 2 ) for t= 1,2. 

Proof. Clearly, the first part of (1) follows from the Propositions 14.31 15.11 and the 
second assertion is a well-known property of the square root of a positive operator 
(cf. Thm. V.3.24 of Ref. [S]). 

Property (2) follows from the inclusion Dom(G 2 ) C Dom(G) = Dom(iV) and 
L e (T)om(N)) C DomfiV 1 / 2 ). □ 

We now apply the sufficient condition for conscrvativity taking as the operator 
G (cf. Ref. [T7]) the self-adjoint operator 

Xu = rp 2 + (pq + qp) + (u) 2 r + 2j)q 2 on Dom(X) = Dom(TV) 

with r > (27) _1 . 
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Proposition 6.2. Suppose that V is twice differentiable and satisfies the growth 
condition \3.'2\) . Then there exists a positive constant b such that 



i 

(6.2) 2Vle(Xu, Gu) + ^ 



V 1/2 L^t 


2 

< b 


X l l 2 u 









for all u G Dom(JV). 

Proof. We first check the above inequality for u <G G£°(R). The vector u clearly 
belongs to the domain of the operators XG, G*X, L* t XLi. Therefore, the left-hand 
side is equal to (it, C(X)u). From (|3.3[) we have then 

(6.3) (u, C(X)u) = -2 (it, ((2 7 r - l)p 2 + u 2 q 2 ) u) 

+ (2rD pp + AD pq + 2{uj 2 r + 2j)D qq ) \\u\\ 2 + i (u, [V(q), X]u) . 

Estimating the commutator as follows (cf. (|3.3p ) 

\i(u,[V(q),X)u}\ = \-(u,(r(pV'(q)+V'(q)p) + 2qV'(q))u))\ 

< 2r\\pu\\ ■ ||V"(g)u||+2||gu|| ■ ||V"(g)u|| 

< {u 1 (r 2 p 2 + q 2 )u) + 2(ii 1 \V'{q)\ 2 u) 

< (u, (r 2 p 2 + (2g 2 v + l)q 2 + 2g 2 v )u) , 

and putting C5 = max { r 2 , 2g 2 r + 1, 2rD pp + AD pq + 2(uj 2 r + 2^)D qq } we find 

\(u,L(X)u)\ < c 5 (u, (p 2 + q 2 + 2) u) 
= c 5 (u, (27V + 3) u) 
< b(u,Xu) . 

as in the proof of Prop. 14.31 Hence, we have now proved the inequality (|6.2p 
for u G C^°(E). The extension to arbitrary u G Dom(A r ) follows by a standard 
approximation argument. □ 

This result yields 

Theorem 1. Suppose that the potential V is twice differentiable and satisfies the 
growth condition Ii3.2\) . Then the minimal semigroup associated with the above 
operators G, Li,Ij2 is Markov. 

Proof. It suffices to apply Theorem 4.4 from Ref. [T7], choosing the positive, self- 
adjoint operator 4> := —Go introduced in the proof of Prop. l^TTl The hypothesis C 
holds by Propositions 16.11 and 16. 2\ and the hypothesis A by Prop. 15.31 Moreover, 
we choose the positive, sclf-adjoint operator G as a sufficiently large multiple of X 
to satisfy 

($5j 1 ,$5, 1 )<(C3a,du) VueDom(X). 

□ 



7. Stationary state 

In order to establish the existence of a steady state of QFP we now start to verify 
the conditions of Theorem VI. 1 in Ref. [T5] , 

Theorem 2. Suppose that 7 > and the potential V is twice differentiable and 
satisfies the growth condition (|3.2p . Then the quantum Markov semigroup (QMS) 
with form generator h6.1\) has a normal invariant state. 
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Proof. We shall apply Theorem IV. 1 from Ref. [T5]- Hypothesis (H) is satisfied 
due to Prop. 15.31 Now choose 

X := rp 2 + (pq + qp) + (uj 2 r + 2 7 ) q 2 , Y := c 6 {p 2 + q 2 ) - c 7 t 

with r > (27) _1 and cq,C7 > 0. Indeed X is clearly positive (cf. Lemma l3.2p . 
Y is bounded below with finite dimensional spectral projections associated with 
intervals ] — 00, A], since it is a translation of a multiple of the number operator. 
In order to check the fundamental inequality 

(7.1) £{X)[u,u] < ~{u,Yu) 

for all u £ Dom(iV) we start from the identity (|6.3p in the proof of Prop. 16.21 and 
estimate the commutator as follows 

\i(u,[V(q),X]u}\ = \-(u,(r(pV'(q) + V , (q)p) + 2qV , (q))u))\ 

< 2r\\pu\\-\\v'(q)u\\+2\\qu\\-\\v'( q )u\\ 

< e (u,(p 2 + q 2 )u) + r —p-(u,\V'( q )\ 2 u) 

< e (u, (p 2 + q 2 )u) + gf—p- (u, (1 + \x\ 2a )u) , 
where we used the elementary inequality 

(i + M 2 ) Q < i + \x\ 2a 

for < a < 1. With the Young inequality we have 

M ~ V [VlXl )S 1/a + 1/(1 -a) ~ aV |X| 
for all 77 > 0. Choosing 77 = (e 2 /r 2 + l) a we obtain 

M*. <^\ x \ 2 + (l-^)(^ 2 + l)" /(1 - a) 
m ~ r 2 + V 1 e 2«/(i-«) 

Therefore we have 

|t («, [V(g), < e (it, (p 2 + q 2 )u) + aeg 2 v (u, q 2 u) 

2 r + l / r 2 + 1 \ /v I 11 ii2 



+ + (i-«) y-^) INI 

< (1 + a<&)e («, (p 2 + q 2 )u) + g 2 v \\u\\ 2 \± e a +a)/{1 L) + ~ T" J ' 

This inequality and (|6.3p give 

(u, < -2 (u, ((2 7 r - \)p 2 + u?q 2 ) u) + (1 + a<#)e («, (p 2 + g 2 )u) 

+ (2rD pp + AD pq + 2(cu 2 r + 2-y)D qq ) \\u\\ 2 

For all r > (2j)^ 1 , we can take an e small enough such that 

c 6 := 2min{ (2jr - 1),oj 2 } - (1 + agl)e > 0. 

Putting 

(1 - a)(r 2 + l) 1 /^-") r 2 + 1 
c 7 := 2rD w + 4Z? P9 + 2(^ 2 r + 2ry)D n + g 2 v "- ^ +a)/(1 _ a) + 9v—^~ 

we find the asserted inequality (|7.1[) : 

{u,C(X)u) < -c 6 (u,(p 2 + q 2 )u)+c 7 \\u\\ 2 
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for u <G C£°(R). The extension to arbitrary u £ Dom(N) follows by a standard 
approximation argument. 

Clearly G is relatively bounded with respect to X by previous results (cf. §4,5). 
Moreover we have (A + X)" 1 (Dom(iV)) = Dom(W 2 ), and hence 

Lt ((A + X)- 1 Dom(iV)) C L f (Dom(iV)) C Dom(JV 1/2 ) = Dom(A 1/2 ) 

for i = 1, 2; A > 1. Therefore Theorem IV.1[T^ can be applied and it yields the 
existence of a normal invariant state. □ 

We remark that the operators X and Y do not commute here, in contrast to 
most examples in Ref. [T5] , 

8. IRREDUCIBILITY AND LARGE TIME BEHAVIOR 

A QMS T on £>(h) is called irreducible if the only subharmonic projections IT in 
h (i.e. projections satisfying 71(H) > H for all t > 0) are the trivial ones or 1. 

If a projection II is subharmonic, the total mass of any normal state a with 
support in p (i.e. such that IL71I = ILt = oTI = a), remains concentrated in II 
during the evolution. Indeed, the state 7l t (o~) at time t then satisfies 

1 = tr (% t (a)) > tr (% t (a)U) = tr (crT t (U)) > tr (an) = tr (a) = 1 . 

As an example, the support projection of a normal stationary state for a QMS is 
subharmonic (cf. Th. II. 1 in Ref. [5D] ). Thus if a QMS is irreducible and has 
a normal invariant state, then its support projection must be 1, i.e. it must be 
faithful. 

In this section we shall prove that the QMS we constructed in fJSJS] is irreducible 
if the strict inequality A > holds. The more delicate limiting case A = is 
postponed to the next section. 

Subharmonic projections are characterized by the following theorem |20j . 

Theorem 3. A projection II is subharmonic for the QMS associated with the oper- 
ators G, Li if and only if its range X is an invariant subspace for all the operators 
Pt of the contraction semigroup generated by G (i.e. Vi > : PfX C X) and 

Li (X n Dom(G)) C X 

for all I 's. 

The application to our model yields the following 

Theorem 4. Suppose that A > 0. Then the QMS T associated with (the closed 
extensions of) the operators G,Li given by i5.1\) and \2. 7[ ) is irreducible. 

Proof. Let II be a subharmonic projection with range X C h. We shall prove that 
either II = or II = 1, and hence T is irreducible. 

The domain of N is Pt invariant because it coincides with Dom(G) which is 
obviously P t -invariant. Since both L\,L,2 map Dom(iV) into Dom(A^ 1 / 2 ) and 
Dom(iV 1 / 2 ) into h, we have 

Li (X n Dom(iV)) C A"nDom(iV 1/2 ) , and 

(8.1) Li (X fl Dom(iV n )) C X n Dom(iV"~ 1/2 ) . 

Then, by the linear independence of L\ and L2 (due to A > 0) 

p(^nDom(JV)) C A" n Dom(iV 1/2 ) , and 

q (X n Dom(AO) C X n Dom(A^ 1/2 ) , 
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and thus, since N = (p 2 + q 2 — l)/2, 

N(XnBom(N)) C A". 

For all n > the resolvent operator R(n; G) = (n — G) _1 maps h in Dom(G) = 
Dom(iV), therefore the operator N n = nNR(n;G) is defined on h. It is also 
bounded, because N is relatively bounded with respect to G and the identity 
GR(n; G)u = —u + nR(n; G)u yields the inequalities 

\\N n u\\ < c\\nGR(n;G)u\\ + c \\nR(n; G)u\\ 

< cn \\u\\ + c(n+ 1) \\nR(n; G)u\\ , 

for all u G h. 

The subspace A is clearly 7V„-invariant because 

/•oo 

R(n; G)= I cxp(-nt)P t dt 
Jo 

maps it into X n Dom(G) = X n Dom(iV) and N(X n Dom(JV)) C A". It follows 
that A is invariant for the operators e~ tNn = Y^k>o(~ tN n ) /k\ (t > 0) of the 
(semi)group generated by N n . 

Notice that, for all u £ Dom(TV) = Dom(G), we have 

\\N n u-Nu\\ = \\N(nR{n;G)u-u)\\ 

< c\\G(nR(n;G)u-u)\\ + c\\nR(n; G)u - u\\ 
= c \\nR(n; G)Gu - Gu\\ + c \\nR(n; G)u - u\\ . 

Thus, by the well-known properties of the Yosida approximations, N n u converges 
towards Nu as n tends to infinity for all u G Dom(iV). It follows then from a 
well-known result in semigroup theory (cf. [21], Chapter IX) that the operators 
e -tN„ conver g e towards the operators e~ tN of the semigroup generated by — N in 
the strong operator topology on h uniformly for t in bounded subsets of [0, +oo[. 
Therefore X is e _tJV -invariant for all t > 0. 

Since the operators e~ tN ,t > arc compact and self-adjoint, it follows that X 
is generated by eigenvectors of N, 

X = Lin { ej | j G J } , 

for some J C No- Here {e.j)j>o is the canonical orthonormal basis of ^ 2 (No). More- 
over, using Xoo := X n n >i Dom(iV"), we have Linje^ | j e J} C X^ C X, since 
ej are the eigenvectors of N. 

By (|8.1j) . Xoo is ii, £2 - hrvariant. Since a and are related to L\,L,2 by an 
invertible linear transformation, X^ is also a, a^-invariant. Therefore, if J is not 
empty (i.e. X ^ { }) taking m = min J we find, iiiimedicitely that Co — & m c m 
belongs to X^ C A". Thus, any efc with fc > belongs to X^ because a tfe e = \f~k\&k 
and A = Aoo coincides with the whole of h. □ 

Corollary 8.1. If A > 0, then all normal invariant states are faithful. 

Proof. The support projection of an invariant state is subharmonic (and non-zero). 
Since the QMS is irreducible, all non-zero subharmonic projections must coincide 
with the identity operator 1. Hence, any normal invariant state must be faithful. 

□ 

Recall the following classical result [22] . 

Theorem 5. Let T be the unital minimal QMS associated with operators G,Lg. 
Suppose that T has a faithful normal invariant state p. Then the vector space of 
fixed points 

T(T) = { x g B(h) | T t (x) = x, Vt > } 
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is an algebra and the invariant state p is unique if and only if !F(T) = CI . 

Remember that the QMS under consideration admits a steady state by Theo- 
rem [5] and by Corollary 18.11 this invariant state is faithful provided that we have 
strict inequality in the Lindblad condition (|1.2p . 

The next result, taken from Proposition 2.3 in Ref. [23] and the Correction in 
Ref. allows us to apply immediately the above theorem and show that the 

QMS converges towards its unique invariant state. Let us introduce first 

N(T) = {xe B(h) | T t (x*x) = T t {x*)T t (x) ,T t (xx*) = T t (x)%(x*)\/t > 0} . 

If T has a faithful normal invariant state, it is easy to show by the Schwarz property 
%(x*)%(x) < T t (x*x), that Af(T) is an algebra. 

Theorem 6. 2.'! .'52 LetT be the unital minimal QMS associated with the operators 
G, Li. Suppose that: i) there exists a domain D c which is a core for both G and G* , 
and ii) for all u £ D c , its image R(n;G)u belongs to Dom(G*) and the sequence 
(nG*R(n; G)u) n >i converges strongly. 
Then 

(1) f(T) = {H,L 1 ,L$,L 2 ,L* 2 }' , 

(2) N(T) C { L\,L\,L 2 ,L\ }' , 

(3) if J-(T) = 7V(T) then for all initial states a, 7l t (o~) converges as t goes to 
infinity towards an invariant state in the trace norm. 

In the above Theorem the set {H, L\, L\, L 2 , LZ,}' denotes the commutant, i.e. 
the set of all operators that commute with H as well as with L\, L\ and L 2 , L\ 
(analogously for {L\, L*, L 2 , L^}') and R denotes the resolvent. Because L\ and 
L 2 are linearly independent, as long as we have strict inequality in (|1.2[) . we see 
that {Li,L 2 }' consists only of operators commuting with q and p. This yields 
W(T) = CI and since CI C T(T) C J\f(T) also CI = T{T) = Af(T). 

To apply Theorem[6]and thus ensure convergence we need to verify the conditions 
(i) and (ii). Condition (i) is obvious from Prop. I5TH (one might pick C™°(IR) as a 
core). Condition (ii) is also easily checked because the operators G, G* and N 
have the same domain by Prop. 15.11 Moreover, their graph norms are equivalent 
(cf. RcmarkEU). Now, for all u G D c = C™(R), since nGR(n; G)u = nR(n; G)Gu, 
the sequence (nR(n; G)u) n >\ converges to u in the graph norm of G. By Rcmark l5.2l 
it is also convergent in the graph norm of G* , i.e. (nG*R(n; G)u) n >\ converges to 
G*u. 

Altogether we have proved the following 

Corollary 8.2. Let 7 > and V £ C 2 (M.) satisfy If A > the QMS 

associated with G and Li has a unique faithful normal invariant state p. Moroever, 
for all normal initial states a , we have 

lim % t (a) = p 

t— >oc 

in the trace norm. 

9. The limiting case A = 

Here we study the case when A = D pp D qq — D pq — 7 2 /4 = and 7 > 0. We start 
with the case V = 0. In this situation we can compare our result to the explicit 
formula for the (unique normalized) steady state in Ref. [5] for zero perturbing 
potential. The kernel of the density matrix of the steady state can be calculated 
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by means of Fourier transform [B] and is given by 
(9.1) 



Poo(x, y) = — 1== cxp (- - 1 [7 2 w 2 (a; + y) 2 + Q(x - yf 
7TV7022 V 47Q22 



. Q12 fx 2 -y 2 
x exp 1 —ilo— — 
Q22 



where the abbreviations are 



Qxi = D pp + u 2 D qq 
Q12 = 2u>'-fD qq 

Q22 = D PP + u 2 D qq + 4:-f(D pq + jD qq ) and 

Q = Q11Q22 — Q12 ■ 

One can see that this becomes a pure state if and only if Q = j 2 oj 2 . We will now 
discuss its implications on the relation between diffusion and friction coefficients. 

Lemma 9.1. Let V = 0, 7 > 0, and the Lindblad condition hold. The steady state 
given by (|9.ip is a pure state, i.e. Q = j 2 ui 2 , if and only if 7 < ui and 

•,2 



(9-2) - D pp D qq -D 2 pq -^/4 

(9.3) D p g = -jD qq 

pp Dqq . 



(9.4) D pp = u?D, 



In this 



case 



D qq = and 
2^/uj z — 7 2 

1 



(9.5) Poc{x,y) = -Vu> 2 -~/ 2 e-( cx +cy > with 

7T 
1 



(9.6) c = -(yV^f+iT). 
Proof. We rewrite 

(9.7) = Q - 7 2 co 2 = (i--L) (Dpp - D qq tu 2 ) 



2 



2 / 2 

(9.8) + JL < D pp + 2D pq — + D qq LU 2 

(9.9) + w(D pp D qq ~D 2 pq -^ 

In the case 7 < u> this directly implies (|9.2[) - (|9.4p , because of 
For 7 = u> the term (|9.8[) is a quadratic polynomial in to. But it has no real zero, 
since D pp D qq - D 2 pq >0 by 
For the case 7 > w we rewrite (|9.7p ~ (|9.9p as 

(9.10) = (D pp + 2D pql + D qq oj 2 f 

(9.11) + 4 7 2 (D pp D qq - D 2 pq - -L 

(9.12) + 7 2 ( 7 2 -^ 2 ), 

which is a contradiction since the summand (|9.11[) is non-negative by (jTT 
The form of D qq and the kernel of the density matrix of the steady state follow by 
straightforward calculations. □ 
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Corollary 9.2. Let V = 0, A = 0, and < 7 < u. Under the conditions (|9.2p , 
(|9.3p and ()9.4|) of Lemma WJ] the semigroup is not irreducible. It admits a steady 
state that is not faithful. 

We now want to interpret this result in terms of the generators. Condition (|9.2p 
implies that the Lindblad operators L\ and L2 are no longer linearly independent 
since Li is zero. The operator L\ has a nontrivial kernel containing functions of 
the type exp(— cx 2 ), with c from (|9.6j) . which are clearly also in the kernel of L\L\. 
Hence, density matrices with kernel (|9.5[) are initially not affected by the dissipative 
part of the evolution. However the Hamiltonian evolution might move them away 
from the kernel. It is exactly conditions (|9.3p and (|9.4p that ensure that H is (up to 
an additive constant) a scalar multiple of L\L\ and thus [H,L*Li] = 0. Of course 
in this case functions of type ()9.5p stay unaffected by the evolution. 

Using the explicit formula (19. ip for the kernel of the density matrix in the case 
V = one can check the non-faithfulness of the steady state also directly: 

Remark 9.3. For p not to be faithful we have to find a u € h with « / 0, such 
that (u,pu) = 0. By positivity of p this is equivalent to pu = 0. Using (|9.ip . such 
a u satisfies 



V V22 V ^ // 

for all x € R. Looking pointwise in x we can drop the factors depending only on x 
and get 



Provided the steady state is not pure, i.e. Q ^ 7 2 w 2 , one can interpret this as 
the Laplace-transform of the term in the second line. Due to the decay of the 
exponential factor it exists for all x £ R and it can not be zero for all x unless u is 
zero. This shows that the steady-state p is faithful except when Q = ~/ 2 lu 2 . In the 
case of a pure steady state the first exponential becomes one and the above integral 
vanishes for all functions that are orthogonal to 



\ 47V22 / 

Note that the complex conjugate of this function is in the kernel of L\, as noted 
earlier. 

The above discussion on the explicit steady state for V = and algebraic com- 
putations of the invariant subspaces for G and L\ showing that they should be 
trivial, lead us to the following 

Conjecture. When A = but V ^ (with V G C 2 {R) satisfying the semi- 

group has a unique faithful normal invariant state and the conclusion of Corollary 
holds. 

The full proof of this conjecture, however, entails a lot of analytical difficulties 
and technicalities on invariant subspaces for strongly continuous semigroups and 
their perturbations and will be postponed to a forthcoming paper. 
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